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Compact QED under scrutiny: it's first order 
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We report new results from our finite size scaling analysis of 4d compact pure U(l) gauge theory with Wilson 
action. Investigating several cumulants of the plaquette energy within the Borgs-Kotecky finite size scaling scheme 
we find strong evidence for a first-order phase transition and present a high precision value for the critical coupling, 
(5t in the thermodynamic limit. 



1. INTRODUCTION 

The nature of the phase transition in com- 
pact QED has been under debate for long time. 
We have addressed this problem in high statis- 
tics runs to reach a final conclusion in that mat- 
ter. An important ingredient of our approach is 
the finite size scaling (FSS) theory a la Borgs- 
Kotecky (BK) first established a long time ago in 
the context of strong first order phase transitions 
[1-3]. According to BK the finite volume parti- 
tion function at temperature (3 in finite volumes 
with periodic boundary conditions (neglecting in- 
tcrfacial contributions) has the remarkably simple 
form Z = e~ vf ^ = e~ vfl ^ + e -vf 3 <J3)+HX) . 
The functions f\{j3) and /2(/3) denote bulk free 
energy densities in the two coexisting phases 1 
and 2. X stands for the asymmetry parameter 
which is nothing but the relative phase weight in 
the probability distribution P{E). 

A heuristic extension of the BK ansatz to weak 
first order transition was demonstrated for the 
3d 3-state Potts model [4] . The conclusion of our 
work is based on a validation of BK by pertur- 
bative analysis as well as independent ab initio 
determinations of the gap characteristics. 

2. SIMULATION DETAILS 

We consider 4d pure U(l) gauge theory with 
Wilson action S = —f3J2n u>^ cos (^iJ.v(n)), where 
(3 represents the Wilson coupling and 9 ay {n) the 
plaquette angle. We use a lattice of volume V = 
L 4 with periodic boundary conditions. 

We have implemented three different algo- 



rithms for generating the U(l) gauge field con- 
figurations: (a) a local Metropolis (Metro), up- 
dating each link separately, (b) a global hybrid 
Monte Carlo algorithm (HMC) and (c) a com- 
bination of the multicanonical and the hybrid 
Monte Carlo algorithm (MHMC). For details we 
refer to Refs. [5,6]. The cumulative number of 
generated configurations at each lattice size L is 
> 5 x 10 6 . We measure the number of tunneling 
events (flips) as control parameter for the mobil- 
ity of the algorithms and the integrated plaquette 
autocorrelation time Tint which controls the sta- 
tistical quality of each single Markov chain. Runs 
differing by coupling, algorithm, HMC parame- 
ters or by weight function are considered as inde- 
pendent. Our simulation parameters are listed in 
Table 1. 

3. CUMULANTS 

Based on the plaquette operator 

E = 6^ E «»(M»)). (!) 

n , v > fj, 

we consider the following cumulants: 
C V (P,L) = 6V ((E 2 ) - (E) 2 ) , 



Ui{fl,L) = 




In addition to their derivatives with respect to (3 
we measure higher derivatives of the free energy 
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Table 1 

Simulation details. The HMC subscript denotes 
length of trajectory. 



L 





algorithm 


#confx 10 b 


# flips 


T int 


(5 


1 


001700 


Metro 


11.20 


21 170 


104(2) 




1 


001500 


Metro 


9.80 


18300 


104(2) 




1 


001600 


HMC 


2.48 


3577 


130(5) 




1 


001772 


MHMC 


4.90 


7444 


102(2) 


jj — 


1 


007370 








304( 16) 




1 


007370 


HMC 


125 


2 274 


1256(196) 




1 


007370 


HMC 2 


1 25 


554 


649(73) 




1 


007370 


HMC 4 


1 25 


698 






1 


007370 


HMC 6 


1 25 


846 


390(33) 




1 


007370 


HMCq 


1.25 


907 


339(27) 




1 


007370 


HMC^q 


1.25 


921 


328(26) 




1 


007370 


HMC^ 


1.25 


929 


345(28) 




1 


007370 


HMC^2 


1.25 


911 


363(30) 




1 


007370 


HMCn 


1.25 


932 


355(29) 




1 


007370 


HMC 14 


1.25 


854 


382(33) 




1 


007370 


hmc 16 


1.25 


856 


399(35) 




1 


007370 


HMC 


1.44 


1058 


379(30) 




1 


007337 


MHMC 


6.36 


5179 


240(7) 


10 


1 


009300 


Metro 


4.37 


1770 


784(52) 




1 


009400 


Metro 


7.44 


3104 


775(35) 




1 


009300 


HMC q 


1.00 


294 


948(144) 




1 


009300 


HMC ii 


1.00 


350 


897(132) 




1 


009300 


HMC 15 


1.00 


353 


1060(170) 






UU^oUU 


HMC17 


1.00 


344 


831(118) 




1 


009300 


HMC iQ 


1.00 


340 


894(132) 




1 


009300 


MHMC 


2.61 


1318 


412(56) 


12 


1 


010143 


Metro 


3.62 


569 


2406(304) 




1 


010143 


Metro 


5.88 


916 


2058(189) 




1 


010143 


Metro 


1.75 


315 


2576(486) 




1 


010143 


Metro 


1.87 


308 


2098(345) 




1 


010143 


MHMC 


1.30 


403 


734(86) 




1 


010143 


MHMC 


2.18 


702 


689(60) 


14 


1 


010598 


Metro 


3.90 


215 


5980(1150) 




1 


010600 


Metro 


6.52 


395 


7480(1240) 




1 


010568 


HMC 


0.57 


26 


12500(9000) 




1 


010568 


MHMC 


0.83 


169 


1070(190) 




1 


010568 


MHMC 


3.80 


725 


1380(130) 


16 


1 


010753 


Metro 


5.42 


75 


22900(9000) 




1 


010800 


Metro 


5.55 


93 


25400(6800) 




1 


010753 


MHMC 


0.60 


66 


1980(560) 




1 


010753 


MHMC 


0.63 


77 


1850(500) 




1 


010753 


MHMC 


1.63 


151 


1770(290) 




1 


010753 


MHMC 


3.39 


424 


1800(200) 


IS 


1 


010900 


MHMC 


0.38 


13 


9900(7900) 




1 


010900 


MHMC 


0.63 


20 


10000(6800) 




1 


010900 


MHMC 


0.98 


58 


5800(2200) 




1 


010900 


MHMC 


3.79 


272 


4700(820) 



density {-l) n+1 n n {{3 , L) = 9 "ffi L) ■ Introducing 
the central moments [i n = V n ^ 1 ((E - (E)) n ) we 
can write them as 

«3 = M3, 

k 6 = /i 6 - 15V>2M4 - 10V/4 + 30V 2 14- 

For each of the ten cumulants the location 
(f3 K , n) of its rightmost extremum is determined 
by reweighting the measured probability distribu- 
tion P(E) to different couplings (3. To calculate 
the estimates of our cumulants at each lattice size 
L we proceed in two steps: i) we determine the 
error of each individual run performing a jack- 
knife error analysis by subdivision of the run into 



ten blocks; ii) we calculate the final result by % 2 - 
fitting these individual results to a constant. 

4. FINITE SIZE SCALING 

Let us consider the expansion of the pseudo- 
critical coupling for C v ; in the BK representation 
it is given by 

(3cv(V) = pcv(oo) + £ B k V- k . (2) 
fc=i 

In order to expose systematic effects in the fit 
parameter 0cv (°°)> we varv both, the fit range 
within L min < L < 18 and the truncation pa- 
rameter k max . Table 2 displays a remarkable 
stability pattern both for pev and B\ support- 
ing the validity of the ^expansion. Averag- 
ing the best fit couplings to a constant we obtain 
0cv(oo) = 1.0111310(62). Performing the same 
analysis for all ten cumulants yields an average 
infinite volume transition coupling 

[3 T = 1.0111331(21). (3) 



Table 2 

Transition couplings Pcv(oo) fitted to Eq (2). 
Best fits are in bold face letters. 



Lmin 


kmax 


Xdof 




Bi 


14 


l 


1.03 


1.0111241(13) 


-18.95(14) 


12 


l 


1.09 


1.0111144(55) 


-18.24(21) 




2 


0.19 


1.0111315(57) 


-19.96(53) 


10 


l 


12.7 


1.0110945(147) 


-17.18(37) 




2 


0.13 


1.0111283(25) 


-19.63(15) 




3 


0.21 


1.0111319(62) 


-20.06(65) 


8 


1 


108 


1.0110474(349) 


-15.33(50) 




2 


2.14 


1.0111159(69) 


-18.70(25) 




3 


0.11 


1.0111309(25) 


-19.94(17) 




4 


0.21 


1.0111309(22) 


-19.94(10) 


6 


1 


970 


1.0109389(913) 


-12.38(56) 




2 


37.1 


1.0110792(218) 


-16.84(41) 




3 


1.25 


1.0111199(55) 


-19.02(22) 




4 


0.10 


1.0111316(11) 


-20.02(6) 



Analogously the expected scaling of the max- 
ima of the specific heat yields a prediction for the 
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infinite volume gap G 



6V 



fe=i 

G = 0.026721(59). 



(4) 



Scaling of the Binder cumulant yields 

oo 

U 4 ,min = U + J2 A kV~ k 
fe=l 

U = -5.816(27) 10~ 4 . (5) 



From Bi — —^p- we can derive an asymmetry 
ln(X) = 3.21(10). 

5. CONSISTENCY CHECKS 

An independent leading order perturbative lat- 
tice calculation confirms the value of hi(X) with- 
out relying on the validity of the BK ansatz [7]. 
In the Coulomb phase the partition function can 
approximately be written as 



2vr 



\(v-i) 



^2(1 - cosp M ) 



leading to the free energy F 2 = In Z 



\(v-i) 



ln(f|) + 21nL - E P ' 1^2(1 - cos^). The 
summation over all momenta can be done 
for asymptotically large L as E P ' m E A1 2(l — 
cosp M ) = aV + 2\nL - b + 0(L~ 2 ) with pa- 
rameters a = 1.999708 and b = 1.701216 that 
can be numerically computed to arbitrary preci- 
sion. With F 2 (0,L) = Vf 2 (0) + AF 2 (P,L) we 
obtain for the Coulomb phase finite size correc- 
tion AF 2 = 6+f ln(27r/ e 2j) = 3.15(8). e R is taken 
from a very accurate measurement of the renor- 
malizcd fine structure constant at the phase tran- 
sition a T = e 2 R /4ir = 0.19(1) [8]. Note that the 
logarithmic correction cancels out for the 4d sym- 
metric system with periodic boundary conditions. 

One can argue that the leading finite size effects 
in the confined phase due to ++ gaugeballs and 
string states are at least 3 respectively 6 orders 
of magnitude smaller than AF 2 on a lattice as 
small as L = 16 [7]. Thus we neglect these con- 
tributions and obtain, in perfect agreement, an 
asymmetry parameter ln(A") = AF 2 = 3.15(8). 



Furthermore our very accurate value of (3t 
(Eq (3)) admits a direct measurement of la- 
tent heat and Binder cumulant, performing 
metastable simulations at ftp on lattice sizes up 
to 32 4 [9] . Denoting the energy peaks of the prob- 
ability distribution in the confined and Coulomb 
phases by E\(L) and E 2 (L) we fit their contin- 
uum limit values Ei = Ei(oo) and find 

G = E 2 -E l = 0.026685(54), (6) 
U = - { 2 „ V = -5.777(16) 10~ 4 . (7) 

Both values are in perfect agreement with the BK 
results from Eqs (4,5). 

6. SUMMARY AND CONCLUSION 

All cumulants investigated in our high statistics 
analysis at L=6, 8, 10, 12, 14, 16, 18 can be described 
by BK first-order FSS. Ab initio measurements 
confirm the FSS results for the infinite volume 
gap G, the Binder cumulant U and the asymme- 
try ln(A). The non vanishing values for G and U 
lead to the conclusion that the phase transition 
in compact 4d £7(1) theory with Wilson action is 
first-order. 
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